o 

(N 



in 

(N 



X 



o 



t 



Mon. Not. R. Astron. Soc. 000, 000-000 (DODO) Printed 26 July 2011 (MN 1^1^ style file v2.2) 

Modelling interaction of relativistic and 
non-relativistic winds in binary system PSR 
B1259-63/SS2883^ - II. Impact of magnetization and 
anisotropy of the pulsar wind 



Bogovalov^ S.V.,Khangulyan^ D.,Koldoba^ A.V.,Ustyugova^ G.V.,Aharonian^'^ F.A. 

^^ ' ^Moscow Engineering Physics Institute (state university), Moscow, Russia, 

■^Institute of Space and Astronautical Science/ J AXA, Sagamihara, Japan, 



O , , '^Keldysh Institute of Applied Mathematics RAS, Moscow, Russia, 

Q ' ^Dublin Institute for Advanced Studies, Dublin, Ireland, 

rr\ ^ Max- Planck- Institut fiir Kernphysik, Heidelberg, Germany. 



cn . ABSTRACT 

00 



^^ ' In this paper, we present a numerical study of the properties of the flow 

produced by the collision of a magnetized anisotropic pulsar wind with its 
environment in binary system. We compare the impact of both the magnetic 
field and the wind anisotropy to the benchmark case of a purely hydrodynam- 
ical (HP) intera c tion o f isotropic winds, which has been studied in detail by 



Bogovalov et al 



(|2008l ). We consider the interaction in axisymmetric approx- 
imation, i.e. the pulsar rotation axis is assumed to be oriented along the line 
between the pulsar and the optical star and the effects related to the pulsar 
orbiting are neglected. The impact of the magnetic field is studied for the 
case of weak magnetization (with magnetization parameter a < 0.1), which 
is consistent with conventional models of pulsar winds. The effects related to 
anisotropy in pulsar winds are modeled assuming that the kinetic energy fiux 
in a non-magnetized pulsar wind is strongly anisotropic, with the minimum at 
the pulsar rotation axis and the maximum in the perpendicular direction. We 
show that, although both considered effects change the shape of the region 
occupied by the terminated pulsar wind, their impact appears to be small. 

Star SS2883 in what follows will be refereed as LS2883 according to the standard catalog of Luminous Stars in the Southern 



Milky Way I Stephenson &: Sanduleak 



197ll) 
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In particular, for the magnetizatfon of the pulsar wind below 0.1, the mag- 
netic field pressure remains well below the plasma pressure in the post-shock 
region. Thus, in the case of interaction of a pulsar with the stellar wind envi- 
ronment (opposite to the case of plerions, i.e. the pulsar interaction with the 
interstellar medium, when the magnetic field becomes dynamically important 
independently on the wind magnetization) the HD approach represents a fea- 
sible approximation for numerical modelling. 

Key words: HD - shock waves - pulsars: binaries. 



1 INTRODUCTION 

Pulsars lose their rotation energy through relativistic winds, the interaction of which with the 
surrounding environment leads to the formation of so-called Pulsar Wind Nebulae (PWNe). 
PWNs are bright in X- and gamma ray regi ons filled with ultrarelativistic electrons acceler - 



ated at the pulsar wind termination shock (JRees fc Gunn 



Kennel fc Coroniti 



1984af ). 



While Crab Nebula is the most famous (although rather at ypical) representative, t he recent 



X-ray and TeV gamma ray observations have shown (see e.g. lGaensler fc Slanell2006l ) that the 
formation of PNWe is a common phenomenon for many other pulsars as well. Importantly, 
the physical conditions, and consequently their observational properties, appear to differ sig- 
nificantly in different PWNs. A very interesting situation arises when a pulsar is located in a 
binary system. In this case the pulsar wind interacts with the wind from the companion star. 
This case, in particular, is realized in the binary system PSR1259-63/LS2883 which consists 
of a ~ 47.8ms pulsar in an elli ptic orbit, with eccentricity e = 0.87, around a massive fast 



rotating 08 optical companion (jjohnston et al. 



1992 



Negueruela et al.ll201ll ). Observations 



show that this system is a strongly variable source of non-therma l emission extending from 
radio to TeV gamma-rays (see e.g. iNeronov fc Chernyakoval 120071 . and references therein). 



e.g. 



'wo other variable 



Paredes et al 



2000 



eV galactic gamma-ray sources, L S 5039 and LS I -|-61 303 (see 



Bosch- Ramon fc Khangulyanll2009l ). are discussed as possible, al 



though less evident cand i dates for being a part of a "binary pulsar" source population 



( Maraschi fc Treves 



19811 : 



Martocchia et al 



I 



2003). LS 5039 consists of an 06.5V star and 



an unidentified compact object in a 3.9 days orbit. This object ha s been detected as a peri- 



odic source of gar nma rays by HESS fiAharonian et al. 



2006 



20051 ) and by Fermi Large Area 



Telescope (LAT) (JAbdo et al.ll2009[ ). In LS I -|-61 303 an unidentified compact source moves 
around a BOVe star in an eccentric orbit with a period of 26.5 days. This source was as well 
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detected in GeV and TeV gamma rays (JAlbert et al. 



2009 



2006 



Abdo et al. 



2009 



Anderhub et al. 



Acciari et al. 



20111). Although, the compact source nature in these objects (black hole 



or neutron star) is not yet firmly established, the binary pulsar scenario provides a possible 
explanation of the broadband non-thermal emission detected from these sources (see e.g. 



Dubus 



20061 ). 



The collision of supersonic winds from two stars located in a binary system results in 
formation of two terminating shock fronts and a tangential discontinuity separating rela- 
tivistic and nonrelativistic parts of the shocked fl ow. It has been found from the detailed 



hydrodynamical modelling ( iBogovalov et al. 



20081 ) that the shocked flow propagates into a 



limited solid angle with rather high velocity. It has been shown that on the scales, which 
characterize the binary sys tem, the shocked pulsar wind material obtains relativistic bulk 



velocities. As proposed by 



Khangulyan et al 



( 120081 ). this should lead to the formation of 



strongly variable non-thermal emission due to the Doppler boosting, thus a detailed study 
of the properties of the post-shock flow is imp ortant for proper interpretation of variable 



X- and gamma rays detected from the system ( Aharonian et £ 



2006 



Uchiyama et al. 



2009; 



Aharonian et al. 



20091 



Abdo et al. 



20n|). 



2005bl : IChernyakova et al. 



In the previous pure HD study by IBogovalov et al.l (120081) several potentially important 



effects have been identifled for the future research. In particular, the development of Kelvin- 
Helmholtz instability at the contact discontinuity and orbital motion of the pulsar should 
enhance the mixing of pulsar and stellar wind materials on a spacial scales slightly exceeding 



the size of the binary system (IBosch-Ramon &: Barkovl 



time of VHE electrons (see e.g. iKhangulyan et al.ll2007l ) this mixing should not signiflcantly 



201ll ). Given relatively short cooling 



affect the variable gamma ray component, which is expected to be produced inside the 
system. On the other hand, several other effects related to impact of magnetic fleld and 
anisotropy of the pulsar wind, which are almost unavoidable features of the pulsar wind, 
appear to be important already on the scales of binary system. The importance of this effects 
can be illustrated by the impact of the magnetic fleld on the interaction of the pulsar wind 
with interstellar medium. Even a weak magnetization of the pulsar wind can result in a 



dramatic modiflcation of the flow in the post 



Bogovalov et al. 



2005 



Bucciantini et al. 



shock region (jKomissarov fc Lyubarsk 



3 



2003 



20051 ). In this paper we present a study of the role 



of the magnetic fleld and anisotropic energy flux on the dynamic of the shocked material in 
binary pulsar systems. 
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2 PROPERTIES OF THE PULSAR AND STELLAR WINDS OF PSR 
B1259-63/SS2883 

The separation distance between the pulsar and optical star in PSR-B1259-63/LS2883 
changes by more than order of magnitude depending on the orbital phase. Consequently, 
the physical conditions at the interaction point can vary significantly. While the processes 
related to particle accele ration and radiation may be strongly affected by change of the sep- 



aration distance (see e.g. iKhangulyan et al.ll2007l ). the regime of hydrodynamical interaction 



of the stars depends weakly on this parameter. Indeed, since the interaction of the winds 
depends on the ratio of the ram pressures of the winds, collision regime remains unaffected 
by the change of the separation distance, unless the stellar wind u ndergoes a significant ac - 



celeration or pulsar wind suffers a strong Compton drag (see e.g. IKhangulyan et al.l 120071 ). 
Finally, the interaction at different orbital phases may differ significantly due to presence of 
dense stellar disk in the system. 

2.1 Parameters of the pulsar wind 

Spin-down luminosity of the pulsar PSR B1259-63 is Lsd = 8 x lO'^^ergs"^. This energy 
is released in the form of electromagnetic and kinetic energy fiuxes carried away by pulsar 
wind. Magnetic field in the wind may be represent as superposition of two components with 
essentially different spatial dependencies: the poloidal magnetic field component decreases 
with radius as r~^, while the toroidal magnetic field has a weaker dependence, namely 
proportional to r~^. Thus, at distances larger than light cylinder radius the toroidal magnetic 
field strongly dominates over the poloidal one. Therefore, only the toroidal magnetic field 
component has been taken into account in the modelling. To keep axial symmetry in the 
model the axis of rotation of the pulsar is assumed to be directed along the line connecting 
the pulsar and optical star. Figured] shows a schematic representation of the wind interaction 
geometry. 

The magnetic field strength is characterized by the magnetization parameter a, which 
is the ratio of the magnetic field energy to kinetic energy fiuxes. Conventional theories of 
plasma generation in pulsar magnetospheres predict a highly magnetized wind with a ^ 1 
at the light cylinder. In this case the energy fiux from pulsar is not isotropic and is expected 
to have the following form 

_ = -,„»c=_ + P„W(9). (1) 
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Here Pnsin^(^) is the Pointing flux obtained in the frameworks of the spht-monopole model 



(JBogovalovl I1999I ) : ^ is the particle flux in the wind; and 6 is the angle between the pulsar 



rotational axis and direction of the flow. Obviously, in this case the ratio of the Poynting 
flux to the kinetic energy density flux varies with 9. It is convenient to define a as the ratio 
-Po/(7o"^c^^) which equals to the ratio of the Poynting flux to the kinetic energy flux in 
the equatorial plane of the pulsar wind. 

Interpretation of observations usually leads to the conclusion that the wind energy flux 
at large distances from the p ulsar is concentrated in the particle kinetic energy, i.e. a <^ 



1 ( Kennel fc Coroniti 



1984bl ). Thus, a transformation of the Poynting flux to the pulsar 



wind bulk motion is to occur in the pulsar wind zone, although we note that the physical 
mechanism responsible for this transformation still remains unknown. Assuming that the 
most of the Poynting flux is transformed into the kinetic energy flux, Eq. ([1]) can be rewritten 



as 



dL r,dN r,,^. r,dN „ r, , ^. ,. 

:^ = ^omc^Z^ + 7max sm2 ^mc^-— + Psm^^ , 2 

ml ml ml 



where we assumed that the particle flux -ttj is isotropic. Then the a para meter equals to the 



ratio P /iimi^^mc^^). We should note that the original definition of a in 



Kennel &: Coroniti 



(Il984bl ) does not differ essentially from the one used here given that 7max ^ 7o- According 
to this definition a parameter is the fraction of the initial Poynting flux remaining in the 
electromagnetic form at the wind termination shock. 

Thus, the requirement of a kinetically dominated pulsar wind leads to the conclusion 
that the wind should be strongly anisotropic and weakly magnetized. To study the impact 
of these effects, below we consider two different models: (i) weakly magnetized isotropic 
wind; and (ii) highly anisotropic unmagnetized wind. 

2.2 Parameters of the stellar wind 

Given the fast rotation of LS2883, the stellar wind is believed to have a strong angular 
dependence. Namely, it is expected that the stellar outflow consists of two parts with com- 
parable mass loss rates: an almost isotropic polar wind and a dense Keplerian disk in the 
plane perpendicular to the star rotation axis. The mass loss rate and velocity profiles of 
stellar winds is a subject of future dedicated studies. We note however that although the 
precise measurements presently are not available, the prop erties of the win d can be esti- 



mated based on the luminosity and temperature of the star (jVink et al.ll2000l ). The physical 
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Figure 1. The scheme of the wind interaction and the post shoclc flow. The flow from optical star is terminated by the 
nonrelativistic shock front (thick dashed line). The flow from the pulsar is terminated by the relativistic shock front (thick 
solid line) . Nonrelativistic and relativistic post shock flows are separated by the contact discontinuity (dashed-dotted line) . The 
symbol © shows the direction of the toroidal magnetic field which has been assumed to exist only in the region occupied by 
pulsar wind. 

characteristics of the optical star L S2883 have been recent ly significantly revised towards 
higher temperature and luminosity (JNegueruela et al.ll201ll ). This leads consequently to a 
higher mass loss rates. In particular, t he mass loss rate of the polar wind is expected to be 



at the level of M = 6 x lO^^Moyr-^ flNe 



meruela et al. 



the rate assumed in the previous study by 



20111). w hich exceeds significantly 



Bogovalov et al.l (120081 ). Thus, the expected ratio 



of the winds ram pressures, r] parameter, for the case of the interaction of the pulsar wind 
with the stellar polar wind is as follows: 



^ -sd ^ 5 ^ iQ-2 -_^ '" (3) 

' McK, V6 X lO-^Moyr-i; VlSSOkms^V ' ^' 

where the wind velocity has been normali zed to the value of the wind termination velocity 
l^nf = 1350 ± 200 km s"^ JmcCoUuJ |i993|) . 



jroad range of variatio n of the r] parameter has been discussed in the 



Since a rather 
previous study by 
r] = 0.1, which are similar to the estimate given by Eq. ([3]). 



Bogovalov et al.l (120081 ). here we adopt two values of 77: r] = 0.05 and 



3 BASIC ASSUMPTIONS. 

The evolution of the fiow due to the change of the star separation distance is assumed 
to be adiabatic, i.e. we consider the collision of the winds as a steady process. All factors 
resulting in time- dependent effects were ignored or intentionally suppressed. In particular. 



the tangential instability arising at the inter 



flows in the post shock region is suppressed ( JBogovalov et al. 



ace between the relativistic and nonrelativistic 



20081). The rotation axis of the 
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pulsar was assumed to be co-directed with the hue connecting the pulsar and the optical 
star, thus the flow in the computational domain was considered as axisymmetric in respect to 
this line. It implies that all effects related to the orbiting of the pul sar are ignored, although 
these effects may appear to be very important on larger scales (JBosch-Ramon fc Barkov 
201lk 

The Lor entz factor of the particles in t he pulsar wind is expected to be very high. 



7n 



10^ (JBogovalov fc Khangoulian 



20021 ). Remarkably, the post shock flow practically 



does not depend on 7max provided that 7max 3> 1. Almost all the energy of the bulk motion 
is transformed into "thermal" energy of plasma, thus it is possible to neglect the rest mass 
energy of the electrons compare with the thermal energy and the particle density is not 
present in the equations describing the motion of the relativistic pulsar wind in the post 
shock region. 



4 METHOD OF SOLUTION 

A special numerical method was used for the integration of the tran sient HP (hydrodynam i- 



cal) and RMHD (relativistic magnetohydrodynamical) equations (jUstyugova et al. 



19991). 



The integration of the HD equations was performed in the region limited by the non- 
relativistic termination shock and the contact discontinuity. The RMHD equations were 
integrated in the region located between the relativistic termination shock and the contact 
discontinuity. The plasma flow outside these two regions corresponds to a supersonic pulsar 
and a stellar wind, i.e. it has been considered as known. The mo delling has been perfor med 
with the use of an adaptive computational mesh (see for details iBogovalov et al.ll2008l ). In 
the region occupied by the shocked relativistic plasma the RMHD equations in cylindrical 
coordinates are 



On'-) 1 drn'^Vr dwyVz 
dt r dr dz 



0. 



(4) 



at + -rT/ '" + '''-'"' + <" + y ' 



+ 



+ — (w + /i^)7V^z 
d^'^jw + h^)v, 1 d . 2^ 2 

h -TT'^iw + h )-f VrV^ + 

at r or 






(5) 



d [ h^ 



0. 



(6) 
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— W{W + /i') - (P + y) J + -^^(^ + h^HVr + 



d 



0. 



dh^ dh-^Vr d'-yVzh 



(7) 
(8) 



dt dr dz 

Here p, w = |p, n and h are the pressure, enthalpy density, particle density and magnetic 
field, respectively (all quantities are in the local rest frame); Vr, and Vz are the components 
of the flow velocity. 

In the region occupied by the shocked nonrelativistic plasma the equations are 



dp 1 drpvr dpVz 
dt r dr dz 







dpVr 1 d / 2 \ 9 

-W^rd-/^P'^^P)^d-zP'' 



p 
r 



dpv 



1 d 



d 



dt +-rd-/^^^' + d-z^'''^^P 



(9) 
(10) 

(11) 



d f pv'^ 
di \~Y' 



1 d ( pv"^ 
r dr ^ \ 2 



P] + 



d I pv^ \ 



(12) 



Here p and e are the densities of the mass and the thermal energy, respectively, and p = 2/3e 
is the pressure. The computational method used in this work has been developed earlier for 
the solution of the equation describing the inter action of a magrietized relat ivistic pulsar 



wind with the interstellar medium (see for details 



Bogovalov et al 



2005 



20081) 



5 RESULTS 

5.1 Impact of magnetic field 

To investigate the impact of the magnetic fleld on the process of the winds collision, a se- 
quence of calculations are performed for different values of the a parameter and a flxed value 
of ?7 = 0.1. The range of considered a parameter has been chosen to be a = 0.003 — 0.1, which 
corresponds to the conventional magnetization of the pulsar wind. As mentioned above, th e 



solution in the hydrodynamical limit (a = 0) was obtained earlier by lBogovalov et al.l (12008[ ). 



In Figure |2] the total pressure distribution (i.e. both gas and magnetic fleld pressures 
p + h'^/2) together with stream lines are shown for different pulsar wind magnetizations. 
A similar distribution obtained for the case of an unmagnetized pulsar wind is shown in 
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Figure 2 . Distribution of total pressure and stream lines in the post shock region for two different values of the magnetization 
of the pulsar wind: a = 0.03(left panel) and cr = 0.1 (right panel). The ratio of the winds ram pressures was assumed to be 
ri = 0.l. 




-0.4 



Figure 3. Distribution of total pressure and stream lines in the post shock region for the unmagnetized pulsar wind, i.e. cr = 0. 
The ratio of the winds ram pressures was assumed to be ?7 = 0.1. 

Fig. [3l A comparison of these figures shows that the impact of the magnetic field on the 
post shock flow is rather weak. The magnetic field coUimates the post shock flow toward the 
axis. This happens due to th e pinching of the flow by the toroidal magnetic field (see e.g. 
Bogovalov fc Tsinganoslll999l ). 



5.2 Anisotropy of the pulsar wind 

As it has been described above, we study the impact of an anisotropic energy flux in the 
pulsar wind under the approximation of an unmagnetized wind. Thus, we have introduced 
a polar angle dependence of the pulsar wind bulk Lorentz factor: 7 = 7o + 7max sin^ 6*, where 
7o is the initial Lorentz factor of the wind in the direction of the rotation axis; and 7max 
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a 0.2 0.4 0.6 O.S r 1 



0.2 0.4 0.6 0.8 



Figure 4. Gas pressure and stream lines in the post shock region for the isotropic pulsar wind with r) = 0.05 (left panel) and 
anisotropic wind with < -q >= 0.05 and a = 100 (right panel). 



is the wind Lorentz factor at the equator. It is expected that the pulsar wind is strongly 
anisotropic, i.e. 7inax ^ 7o- The level of anisotropy can be characterized by the ratio of 
polar and equatorial Lorentz factors: a = 7max/7o • In the calculations we adopted a value of 
a = 100. The definition of the rj parameter should be modified in the case of an anisotropic 
pulsar wind. Namely, it is natural to assign it to the ram pressure ratio averaged over the 
polar angle, i.e. < i] >. 

In Figure IH a comparison of the pressures and stream lines is shown for two cases: 
isotropic pulsar wind with 77 = 0.05 and anisotropic pulsar wind with < t] >= 0.05 and 
a = 100. It might be seen from the figure that the relativistic terminal shock is located closer 
to the pulsar at the axis. This happens because the momentum flux in this direction is much 
smaller than the averaged one. It might be seen as well, that in the case of an anisotropic 
pulsar wind, the tangential instability is somewhat more pronounced. In particular, this is 
visible on the shape of the pulsar wind termination shock. Beside these minor changes the 
geometry of the flow, in the region down stream farther from the axis, is similar to the flow 
produced by an isotropic winds. 



6 DISCUSSION AND CONCLUSIONS 

The present study of the interaction of a pulsar and a stellar winds within the framework of 
an MHD approach shows that the magnetic field has a rather small impact on the post shock 
flow. The only apparent effect revealed by MHD consideration, compared to a pure HD treat- 
ment, is weak coUimation of the relativistic flow towards the axis. This conclusion is in strong 
contrast with the case of plerions, i.e. when the pulsar wind interacts with the interstellar 



medium ( lBucciantinill2002l : lKhangoulian fc Bogovalovll2003 ; 



Komissarov fc Lyubarsky 



2003 : 
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Bogovalov et al.l 120051 ) . In the case of plerions even a weak seed toroidal field results in ef- 
ficient collimation of the post shock flow. This difference has a rather fundamental reason 
illustrated in Fig. [6l where the distribution of the bulk motion Lorentz factor is shown. As 
it is seen from this figure, the bulk Lorentz factor is increasing in the downstream region, 
i.e. physically the wind collision geometry acts as a nuzzle, providing an efficient bulk ac- 
celeration of the plasma in the post shock flow. The shock pulsar wind velocity fast reaches 
relativistic values downstream the shock. It implies an in efficient collimation of the plas ma, 



3 



which is a rather general effect for relativistic outflows ( iBogovalov &: TsinganosI 119991 1. In 
contrast to this case, if the pulsar wind collides with homogeneous medium (i.e. interstellar 
gas), a flow deceleration occurs in the post sh ock region. This results into a few effects mak- 



ing magnetic collimation to be very efficient (JKhangoulian fc Bogovalovll2003l ). The inertia 



of the plasma decreases with a decrease of the velocity, thus it is easier to change the direc- 
tion of the flow even by a small force. Additionally, since in this case the flow is subsonic, 
the toroidal magnetic fleld increases with the distance from the shock. For this reason even 
a weak initial fleld is amplifled quickly and can provide an efficient collimation of the flow 
towards the axis. 

The impact of the pulsar wind magnetization and anisotropy is investigated for case 
of the collision of the pulsar wind with the non-relativistic wind from a massive star. It 
is shown, that the magnetic fleld leads to weak collimation of the post shock flow toward 
the axis. The effect remains very weak for a broad parameter range, in particular for the 
strongest considered magnetization of a = 0.1. This weak impact of the magnetic fleld on 
the wind collision is the result of the acceleration of the plasma bulk motion due to adiabatic 
cooling. This implies that the calculation of the radiation produced at the wind collision can 
be performed in the hydrodynamical limit in which the toroidal magnetic fleld is calculated 
as a passive scalar not affecting dynamics of the plasma. 

The influence of anisotropy of the pulsar wind appeared as well to be small. Whereas 
its impact can be seen on the shape of the pulsar wind termination shock and the enhanced 
potential rate of instabilities, the general structure of the flow remains almost identical to 
the case of isotropic pulsar wind. Although, we note that in the case of a different orientation 
of the pulsar equatorial plane, a stronger effect can arise. 



12 Bogovalov et. al. 




^^11 



U U.2 (J.4 0.6 0.8 r 1 



0.2 0.4 0.6 0.8 



1 



Figure 5. Lorcntz factor in the post shock region of the isotropic pulsar wind with a = 0.03 (left panel) and a = 0.1(right 
panel) . 
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